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Abstract 

Let B^'^ = {B^'^{t), t £ R_|_} be a bifractional Brownian motion in R'^. We prove 
that B^'^ is strongly locally nondeterministic. Applying this property and a stochastic 
integral representation of B^'^ , we establish Chung's law of the iterated logarithm for 
B^'^ , as well as sharp Holder conditions and tail probability estimates for the local times 
of 

We also consider the existence and the regularity of the local times of multiparameter 
bifractional Brownian motion B^-^ = {B^'^{t), t £ in M.'^ using Wiener-Ito chaos 
expansion. 
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1 Introduction 

In recent years, there has been of considerable interest in studying fractional Brownian motion 
due to its applications in various scientific areas including telecommunications, turbulence, 
image processing and finance. On the other hand, many authors have proposed to use more 
general self-similar Gaussian processes and random fields as stochastic models; see e.g. Addie 
etal. (1999), Anhet al. (1999), Benassi et al. (2000), Mannersalo and Norros (2002), Bonami 
and Estrade (2003), Cheridito (2004), Benson et al. (2006). Such applications have raised 
many interesting theoretical questions about self-similar Gaussian processes and fields in 
general. However, contrast to the extensive studies on fractional Brownian motion, there has 
been little systematic investigation on other self-similar Gaussian processes. The main reasons 
for this, in our opinion, are the complexity of dependence structures and the non-availability 
of convenient stochastic integral representations for self-similar Gaussian processes which do 
not have stationary increments. 

The objective of this paper is to fill this gap by developing systematic ways to study sample 
path properties of self-similar Gaussian processes. Our main tools are the Lamperti transfor- 
mation [which provides a powerful connection between self-similar processes and stationary 
processes; see Lamperti (1962)] and the strong local nondeterminism of Gaussian processes 
[see Xiao (2005)]. In particular, for any self-similar Gaussian process X = {X{t),t € M}, the 
Lamperti transformation leads to a stochastic integral representation for X. We will show 
the usefulness of such a representation in studying sample path properties of X. 

For concreteness, we only consider a rather special class of self-similar Gaussian processes, 
namely, the bifractional Brownian motions introduced by Houdre and Villa (2003), to illus- 
trate our methods. Given constants H S (0,1) and K € (0,1], the bifractional Brownian 
motion (bi-fBm, in short) in R is a centered Gaussian process Bq'^ = {B^'^ {t),t E M+} 
with covariance function 

(1.1) 

and 5^'''(0) = 0. 

Let B^'^ , . . . , bI^'^ be independent copies of Bq'^'^. We define the Gaussian process 
B'^'^ = {B"^^{t), t G M+} with values in M'^ by 

B^^^it) = (<'^(t), . . . , B^'^'it)) , yt G M+. (1.2) 

By (jl.ip one can verify easily that B^'^ is a self-similar process with index HK, that is, for 
every constant a > 0, 

{s^-^(at), t G M+} = {a^^-B^'^(t), t G M+}, (1.3) 

where X = Y means the two processes have the same finite dimensional distributions. Note 
that, when K = 1, B^'^ is the ordinary fractional Brownian motion in However, if 
K ^ 1, B^'^ does not have stationary increments. In fact, fractional Brownian motion is the 
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only Gaussian self-similar process with stationary increments [see Samorodnitsky and Taqqu 

(1994) ]. 

Russo and Tudor (2006) have established some properties on the strong variations, local 
times and stochastic calculus of real-valued bifractional Brownian motion. An interesting 
property that deserves to be recalled is the fact that, when HK = ^, the quadratic variation 
of this process on [0, t] is equal to a constant times t. This is really remarkable since as far as 
we know this is the only Gaussian self-similar process with this quadratic variation besides 
Brownian motion. Taking into account this property, it is natural to ask if the bifractional 
Brownian motion B^'^ with KH = ^ shares other properties with Brownian motion (from 
the sample path regularity point of view). As it can be seen from the rest of the paper, the 
answer is often positive: for example, the bi-fBm with HK = ^ and Brownian motion satisfy 
the same forms of Chung's laws of the iterated logarithm and the Holder conditions for their 
local times. 

The rest of this paper is organized as follows. In Section 2 we apply the Lamperti 
transformation to prove the strong local nondeterminism of B^'^ . This property will play 
essential roles in proving most of our results. In Section 3 we derive small ball probability 
estimates and a stochastic integral representation for B^'^ . Applying these results, we prove 
a Chung's law of the iterated logarithm for bifractional Brownian motion. 

Section 4 is devoted to the study of local times of one-parameter bifractional Brownian 
motion and the corresponding A^-parameter fields. In general, there are mainly two methods 
in studying local times of Gaussian processes: the Fourier analysis approach introduced by 
Bcrman and the Malliavin calculus approach. It is known that, the Fourier analysis approach 
combined with various properties of local nondeterminism yields strong regularity properties 
such as the joint continuity and sharp Holder conditions for the local times [see Berman 
(1973), Pitt (1978), Geman and Horowitz (1980), Xiao (1997, 2005)]; while the Malhavin 
calculus approach requires less conditions on the process and establishes regularity of the 
local times in the sense of Sobolev-Watanabe spaces [see Watanabe (1984), Imkeller et al. 

(1995) , Eddahbi et al. (2005)]. In this paper we make use of both approaches to obtain more 
comprehensive results on local times of bifractional Brownian motion and fields. 

Throughout this paper, an unspecified positive and finite constant will be denoted by 
c, which may not be the same in each occurrence. More specific constants in Section i are 
numbered as c. ^ , c. 2 , • • ■ ■ 

2 Strong local nondeterminism 

The following proposition is essential in this paper. From its proof, we see that the same 
conclusion holds for quite general self-similar Gaussian processes. 

Proposition 2.1 For all constants < a < b, Bq'^ is strongly locally ip-nondeterministic 
on I = [a,b] with ip{r) = r^^^ . That is, there exist positive constants c^^ and ro such that 
for all t E I and all < r < min{t, ro}, 

Vav(^B^'^{t)\B^'^{s) ■.seI,r<\s-t\<ro)> c,^, ip{r). (2.1) 
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Proof We consider the centered stationary Gaussian process Yq = {YQ{t),t £ M} defined 
through the Lamperti's transformation [Lamperti (1962)]: 

= e-^^*5^'^(e*), VteM. (2.2) 

The covariance function r{t) := E{YoiO)Yo{t)) is given by 

<2HK 



- A. pHKt 
2^ 



+ 1)" -|e*-l 



2Ht\K H e"*P'^^ 



(2.3) 



Hence r{t) is an even function and, by ()2.3p and the Taylor expansion, we verify that r{t) = 
C»(e-^*) as t ^ oo, where (3 = min{iJ(2 - K),HK]. It follows that r(-) G L'^{M). On the 
other hand, by using (|2.3p and the Taylor expansion again, we also have 

r{t)^l-^\t\^^^ as t^O. (2.4) 

The stationary Gaussian process Yq is sometimes called the Ornstein-Uhlenbeck process 
associated with Bq'^ [Note that it does not coincide with the solution of the fractional 
Langevin equation, see Cheridito et al. (2003) for a proof in the case K = 1]. By Bochner's 
theorem, Yq has the following stochastic integral representation 

YQ{t)= [ e'^^W{dX), VtGM, (2.5) 



where is a complex Gaussian measure with control measure A whose Fourier transform is 
r(-). The measure A is called the spectral measure of Y. 

Since r(-) € L^(R), so the spectral measure A of y has a continuous density function 
/(A) which can be represented as the inverse Fourier transform of r(-): 

/(A) = - / r{t)cos{tX)dt. (2.6) 
^ Jo 

It follows from ()2.4p . (12. 6p and the Tauberian theorem due to Pitman (1968, Theorem 5) [cf. 
Bingham et al. (1987)] that 

/(A) ~ c,^, |^|-(i+2//A') as A ^ oo, (2.7) 

where 2 > is an explicit constant depending only on HK. Hence, by a result of Cuzick and 
DuPreez (1982, Lemma 1) [see also Xiao (2005) for more general results], Iq = {^o(*))* £ ^} 
is strongly locally (/7-nondeterministic on any interval J = [— T, T] with (p{r) = r"^^^ in 
the sense that there exist positive constants 5 and 3 such that for all t € [—T, T] and all 
r G (0, |t| A5), 

Yar{Yo{t)\Yo{s) : s e J, r < \s - t\ < 6) > 3 ip{r). (2.8) 
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Now we prove the strong local nondeterminism of Bq ' on /. To this end, note that 

B^'^it) = t^^Yo{logt) for ah t > 0. We choose vq = a6. Then for ah s,t e I with 
r < \s — t\ < tq we have 

^ < |logs-logt| < S. (2.9) 
Hence it follows from (|2.8p and (|2.9p that for all t G [a, b] and r < tq, 



Ysii\^BQ'^{t)\BQ'^{s) : s G /, r < |s - t| < tq^ 
= Var(^t^^yo(logt)|s'^^>o(logs) :s£l, r<\s-t\<ro^ 

> t^^^Yav(Yo{logt)\Yo{logs) : s e I, r < \s - t\ < r^j (2-10) 

> a^^^Var(^yo(logt)|lo(logs) : s G /, r/6 < | log s - log t| < 

This proves Proposition 12. li □ 

For use in next section, we list two properties of the spectral density /(A) of Y. They 
follow from (j2.7p or, more generally, from (|2.4p and the truncation inequalities in Loeve (1977, 
p. 209); see also Monrad and Rootzen (1995). 

Lemma 2.2 There exist positive constants g and ^ such that for u > 1, 

j AV(A)dA<c,,^2(i-^^) (2.11) 

J |A|<u 

and 

/(A)dA<c,,„^x-2^^. (2.12) 



/ 

J\\ 



'\X\>u 

We will also need the following lemma from Houdre and Villa (2003). 
Lemma 2.3 There exist positive constants c^ ^ and g such that for all s,t G M+, we have 

2i 



B^'^it)-B^'^is)y <c,,,\t-sr^. (2.13) 



3 Chung's law of the iterated logarithm 

As applications of small ball probability estimates, Monrad and Rootzen (1995), Xiao (1997) 
and Li and Shao (2001) established Chung-type laws of the iterated logarithm for fractional 
Brownian motion and other strongly locally nondeterministic Gaussian processes with station- 
ary increments. However, there have been no results on Chung's LIL for self-similar Gaussian 
processes that do not have stationary increments [Recall that the class of self-similar Gaussian 
processes is large and fBm is the only such process with stationary increments]. 
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In this section, we prove the following Chung's law of the iterated logarithm for bifrac- 
tional Brownian motion in R. It will be clear that our argument is applicable to a large class 
of self-similar Gaussian processes. 

Theorem 3.1 Let Bq'^ = {BQ'^{t), t € be a bifractional Brownian motion in R. 

Then there exists a positive and finite constant Cg ^ such that 

. maxtg[o,^] \Bq {t)\ 

imimi ; — - — r^rTTj^ = c, , a.s. o.l) 

^_.o r^^/(loglog(l/r))^^ ^ ' 

In order to prove Theorem I3.H we need several preliminary results. Lemma 13.21 gives 
estimates on the small ball probability of ' . 

Lemma 3.2 There exist positive constants Cg ^ and Cg g such that for all to € [0, 1] and 
2;G(0,1), 

( - ^17^^) ^ ^ {/^o"" ^^^^""^'^ - < < exp ( - . (3.2) 



Proof By Proposition 12.11 and Lemma [2.31 s^e that Bq ' satisfies Conditions (CI) and 
(C2) in Xiao (2005). Hence this lemma follows from Theorem 3.1 in Xiao (2005). □ 

Proposition 13.31 provides a zero-one law for ergodic self-similar processes, which comple- 
ments the results of Takashima (1989). In order to state it, we need to recall some definitions. 

Let X = {X(t), t G R} be a separable, self-similar process with index k. For any constant 
a > 0, the scaling transformation S^^^a of X is defined by 

(S«,a^)(i) = a-^X{at), Vt G R. (3.3) 

Note that X is K-self-similar is equivalent to saying that for every a > 0, the process 
{{Sf^^aX){t),t G R} has the same finite dimensional distributions as those of X. That is, 
for a K-self-similar process X, a scaling transformation Si^^a preserves the distribution of 
X, and so the notion of ergodicity and mixing of S^^a can be defined in the usual way, cf. 
Cornfeld et al. (1982). Following Takashima (1989), we say that a K-self-similar process 
X = {X(t), t G R} is ergodic (or strong mixing) if for every a > 0, a 7^ 1, the scaling transfor- 
mation Sk^q is ergodic (or strong mixing, respectively). This, in turn, is equivalent to saying 
that the shift transformations for the corresponding stationary process Y = {Y{t),t G R} 
defined by Y{t) = e~'^^X{e^) are ergodic (or strong mixing, respectively). 

Proposition 3.3 Let X = {X{t),t G R} he a separable, self-similar process with index k. 
We assume that X{0) =0 and X is ergodic. Then for any increasing function ijj : R+ — > R+, 
we have ¥{E^^^) = or 1, where 

Ek.,iP = sio '■ there exists 6 > such that sup \X{s)\ > t'^ip{t) for allQ <t < 5> . (3.4) 
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Proof We will prove that for every a > 0, the event -E^,?/' is invariant with respective to 
the transformation 8^,0- Then the conclusion follows from the ergodicity of X. 

Fix a constant a > and a 7^ 1. We consider two cases: (i) o > 1 and (ii) a < 1. In the 
first case, since ip is increasing, we have ip{au) > '<p{u) for all u > 0. Assume that a.s. there 
is a > such that 

sup \X{s)\ > e^{t) for all < t < 5, (3.5) 

0<s<t 

then 

snv \a-'' X{as)\=a-'^ sup \X{s)\>t''^{t) for all < t < (3.6) 

0<s<t 0<s<at 

This implies that E^^^ C S~\(yEf^^^Y By the self-similarity of X, these two events have 
the same probability, it follows that P{-E'k^^A5~J(£'k^^)} = 0. This proves that -B^.i/j is 
5K_a-iiivariant and, hence, has probability or 1. 

In case (ii), we have ip{au) < ip{u) for all u > and the proof is similar to the above. If 
Sn.aX € E^K,?/i, then we have X € E'k,^. This implies S~l^[Ef^^^) C -Ek,!/; and again E^^^^ is 
^K^a-invariant. This finishes the proof. □ 



sian processes, we see that Bq ' is mixing. Hence we have the following corollary of Propo- 



By a result of Manuyama (1949) on ergodicity and mixing properties of stationary Gaus- 
n procef 
sition 13. 3[ 

Corollary 3.4 There exists a constant Cg 4 € [0, 00] such that 

(loglogl/t)^ \j,H,Ki .\ i^^. 

Proof We take ^c(*) = c(loglogl/t) and define c^^ = supjc > : PjSK.V'c} = l}- 
It can be verified that (jS.Th follows from Proposition 13.31 □ 

It follows from Corollarv 13.41 that Theorem 13 . 1 1 will be established if we show Cg ^ G (0, 00). 
This is where Lemma 13.21 and the following lemma from Talagrand (1995) are needed. 

Lemma 3.5 Let X = {X{t),t G M} he a centered Gaussian process in M and let S be a 
closed set equipped with the canonical metric defined by 

d{s,t) = \e{X{s) - X{t)fy^^ . 
Then there exists a positive constants ^ such that for all u > 0, 

<expf-^y (3.8) 



4 sup \Xis)-X{t)\>c,Ju+ r .J\ogNd{S,e)de 

[s, teS ' \ Jo 



where Nd{S, e) denotes the smallest number of open d-balls of radius e needed to cover S and 
where D = sup{(i(s,f) : s, t G S*} zs the diameter of S. 

Now we proceed to prove Theorem 13.11 
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Proof of Theorem l3.1l We prove the lower bound first. For any integer n > 1, let r„ = e 

Let < 7 < Cg 3 be a constant and consider the event 

Ar. = [ max \B^'''{s)\ < 7^^rf ^/(log log l/r„)^ 

Then the self-similarity of B^'^ and Lemma 13.21 imply that 

P{An} < exp - ^ log = n'^^s.s/T. (3.9) 

Since Yln=i P{^n} < 00, the Borel-Cantelli lemma implies 

max^gFQ^ 1 \bI^'^{s)\ 
™^ r„^^/(logrog(l/r„))^^ - 
It follows from (j3.10p and a standard monotonicity argument that 

. maxtg[o,,,]^o' {t)\ 
r^Q r^'f^/(loglog(l/r))^^ ^'^ 

The upper bound is a little more difficult to prove due to the dependence structure of 
Bq'^. In order to create independence, we will make use of the following stochastic integral 
representation of Bq'^: for every t > 0, 



= f^^ / e*^^°§*PF(dA). (3.12) 

[ repres 

For every integer n > 1, we take 



This follows from the spectral representation (|2.5p of Y and its connection with Bq'^ 



= n~" and dn = n^, (3.13) 

where /? > is a constant whose value will be determined later. It is sufficient to prove that 
there exists a finite constant Cg ^ such that 

maXggrg t 1 Ldq (sj 

^r^^cS^ ^^^/(iogrog(iA.))^^ - '^^^ ^'-'^^ 

Let us define two Gaussian processes Xn and X„ by 

Xn{t)=t^^ [ e'^^°^'W{dX) (3.15) 



and 



Xn{t) = t"^ j e^^^^'^^WidX), (3.16) 

J|A|^{d„_l,d„l 



respectively. Clearly Bq'^ {t) = Xn{t) + Xn{t) for all t > 0. It is important to note that the 
Gaussian processes Xn{n = 1,2,...) are independent and, moreover, for every n > 1, Xn 
and Xn are independent as well. 

Denote h{r) = r^^ (loglogl/r) . We make the following two claims: 
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(i). There is a constant 7 > such that 



(ii). For every e > 0, 



Vf| max \XJs)\ <j^^ h(tn)\ = 00. (3.17) 
Vp| max |X„(s)| > e/i(t„)l < 00. (3.18) 

^1 U6[0,in]' J 



Since the events in ()3.17p are independent, we see that ()3.14p follows from ()3.17p . ()3.18p and 
a standard Borel-Cantelli argument. 

It remains to verify the claims (i) and (ii) above. By Lemma [3. 2 1 and Anderson's inequality 
[see Anderson (1955)], we have 



max 



Xn{s)\ < l''''h{tn) > P max |i?o^''^(s)| < j''''h{tn) 



se[o,t„] J 1^ sg[o,t 

> exp ^ ^log(nlogn)^ (3.19) 

= (nlogn) "^3,2 /'T^ 

Hence (i) holds for 7 > Cg 2 . 

In order to prove (ii), we divide [0,t„] into pn + 1 non-overlapping subintervals Jnj = 
[anj-i, CLnj], (i = 0, 1, . . . ,pn) and then apply Lemma [331 to Xn on each of Jnj- Let (3 > 
be the constant in (j3.13p and we take J„_o = [0, tnn~^]. After J^j has been defined, we 
take flnj+i = anj(l + n~^). It can be verified that the number of such subintervals of [0, t„] 
satisfies the following bound: 

Pn + 1 < cn^ log n. (3.20) 
Moreover, for every j > 1, if s, t G Jnj and s < t, then we have t/s — 1 < n^^ and this yields 

t-s<sn~^ and log (-) < n"^. (3.21) 



Lemma 12.31 implies that the canonical metric d for the process Xn satisfies 

d{s,t) <c\s -t\"^ foralls,t>0 (3.22) 

and d{0, s) < ct^^ j^-PHK every s & Jno- It follows that Dq := sup{(i(s, t); s, t G Jn 0} < 
ct^^n-P^^ and 

t 

Na{Jn,o,e)<c^j^. (3.23) 

Some simple calculation yields 

/ log Nd{Jn,o, e) de < I y log (^-^^^ j de 



V VM 



''3,8 ^ 



It follows from Lemma 13.51 and (j3.24p that 



max > eh{tn) \ < exp 

SGJn.O 



n 



2I3HK 



2HK / * (3.25) 

log(nlogn)) 

For every 1 < j < Pn, we estimate the d-diameter of Jn,j- It follows from (13.160 that for 
any s, i € Jn,j with s < i, 



E 



(Xn{s)-Xn{t) 



|A|<d„_i 
+ 



^HK giAlogt _ gHK giAlogs 



fWdX 

^ f{\)d\ 



'\\\>dn 

The second term is easy to estimate: for all s,t G J„j, 



(3.26) 



J2<4e^ / /(A)dA<C3,tr^n-2/^^^, 

J\X\>dn 



(3.27) 



where the last inequality follows from (j2.12p . 

For the first term, we use the elementary inequality 1 — cosx < to derive that for all 
s, t G Jn,j with s < t, 

'^^HK _ ^HK^ + 2t^^ .^^(1 - COS (A log -) 



< s 



|A|<d„_i 
2HK 



s ^ Jr Vs/ Jix\<d„ 



f{X)d\ 
X^f{X)dX 



(3.28) 



< .2HK ~2(3HK 

— "^3,10 ''n ) 

where, in deriving the last inequality, we have used ()3.2ip and (j2.1ip . respectively. 
It follows from (I3.26p . (I3.27P and (j3.28p that the d-diameter of Jnj satisfies 

— '-3,11 ''^ ^ 

Hence, similar to ()3.25p . we use Lemma [331 and p.29p to derive 

^2(3 HK 



max X„(s) > eh{tn) } < exp — c 



n 



log(nlogn)) 



2HK / • 



(3.29) 
(3.30) 



By combining ()3.20p . (j3.25p and (j3.30p we derive that for every e > 0, 

00 00 Pn 

Vp| max \Xn{s)\ >e/i(t„)| < VVpI max |X„(s)| >eh{tn)} 
I- se[o,t„] J ^ ^ L seJ„,j J 

n=l n=l 7=0 



00 / 
n=l ^ 



n 



2(3HK 



log (n log n)) 



2HK 



(3.31) 



< 00. 



This proves (13.18P and hence the theorem. 



□ 
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Remark 3.6 Let to € [0, 1] be fixed and we consider the process X = {X{t)^t € M+} defined 
by X{t) = B^'^it + to) — ^^'^{tQ). By applying Lemma [3.21 and modifying the proof of 
Theorem 13. 11 one can show that 

S,i2^^™int ^Hi^/(loglog(l/r))^^ -"^'^^ ^-^--^^^ 

where Cg > 1 is a constant depending on HK only. 

Corresponding to Lemma 13.21 we can also consider the small ball probability of Bq ' 
under the Holder-type norm. For a € (0, 1) and any function y G Co([0, 1]), we consider the 
a-H61der norm of y defined by, 

II I, \y{s) - y{t)\ 

WvL = sup — — — . (3.33) 

s,te[o,i],s^t \s - tr 

The following proposition extends the results of Stolz (1996) and Theorem 2.1 of Kuelbs, 
Li and Shao (1995) to bifractional Brownian motion. 

Proposition 3.7 Let B^'^ he a bifractional Brownian motion in M and a G (0, HK). There 
exist positive constants Cg jg and Cg such that for all e G (0, 1), 

exp ( - cg,,g < p{||5f'^|U < e} < exp ( - Cg,,, (3.34) 

Proof It follows from Theorem 3.4 in Xiao (2005). □ 



4 Local times of bifractional Brownian motion 

This section is devoted to the study of the local times of the bi-fBm both in the one-parameter 
and multi-parameter cases. As we pointed out in the Introduction there are essentially two 
ways to prove the existence and regularity properties of local times for Gaussian processes: 
the first is related to the Fourier analysis and the local nondeterminism property; the second is 
based on the Malliavin calculus and Wiener-Ito chaos expansion. We will apply the Fourier 
analysis approach for the one-parameter case and the Malliavin calculus approach for the 
multiparameter case. 

4.1 The one-parameter case 

Let B^'^ = {B^'^ {t),t G be a bifractional Brownian motion with indices H and K in 
W^. For any closed interval / C M+ and for any x £ W^, the local time L{x,I) of B^'^ is 
defined as the density of the occupation measure fij defined by 

fii{A) = J^lA{B"'^{s))ds, AeB{R'^). 
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It can be shown [cf. Geman and Horowitz (1980) Theorem 6.4] that the following occupation 
density formula hods: for every Borel function g{t, x) > on / x M'^, 



[ g{t,B"'^{t)) dt= [ [ g{t,x)L{x,dt)dx. (4.1) 



Lemma 12.31 and Theorem 21.9 in Geman and Horowitz (1980) imply that if \/{HK) > d 
then B^'^ has a local time L{x,t) := L{x, [0,t]), where {x,t) € R*^ x [0, oo). In fact, more 
regularity properties of L{x,t) can be derived from Theorem 3.14 in Xiao (2005) which we 
summarize in the following theorem. Besides interest in their own right, such results are also 
useful in studying the fractal properties of the sample paths of B^'^ . 

Theorem 4.1 Let B^'^ = {B^'^ {t),t E M} be a bifractional Brownian motion with indices 
H and K in M^. Ifl/{HK) > d, then the following properties hold: 

(i) B^'^ has a local time L{x,t) that is jointly continuous in {x,t) almost surely. 

(a) [Local Holder condition] For every B € ;B(M), let L*{B) = sup^,g]gd L(x, i?) he the 
maximum local time. Then there exists a positive constant c^ ^ such that for all to € 

j.^^^pL%B^^^^^ a... (4.2) 

Here and in the sequel, B{t, r) = {t — r,t + r) and ipi{r) = r^~^-'^'^(loglog 1/r)^^'^. 

(Hi) [Uniform Holder condition] For every finite interval I QR, there exists a positive finite 
constant ^ such that 

V L*{B{to,r)) 

iimsup sup <C42 (4-3) 

r^O toe/ 'P'iv) 

where Lp2{r) = r^-HKd^Yogl/r)"^'^ . 

Proof By Proposition 12.11 and Lemma 12.31 we see that the conditions of Theorem 3.14 in 
Xiao (2005) are satisfied. Hence the results follow. □ 

The following states that the local Holder condition for the maximum local time is sharp. 

Remark 4.2 By the definition of local times, we have that for every interval Q C 

\Q\= [ L{x,Q)dx < L*{Q) ■ ( inax\B^''^{s) - B^'^{t)\Y. (4.4) 

Jbh.K(^q^ \s,teQ ) 

By taking Q = B{tQ,r) in (j4.4p and using ()3.32p in Remark |3.6^ we derive the lower bound 
in the following 

C43 < limsup (-^(^0' < , a.s., (4.5) 
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where C43 > is a constant independent of to and the upper bound is given by (j4.2p . A 
similar lower bound for (j4.3p could also be established by using ()4.4p . if one proves that for 
every interval / C 

liminf inf max 777— , , — < c. ^ a.s. (4.6) 

Theorem 14. II can be applied to determine the Hausdorff dimension and Hausdorff measure 
of the level set Z^ = {t eM+ : B^'^{t) = x}, where x G M"'. See Berman (1972), Monrad and 
Pitt (1987) and Xiao (1997, 2005). In the following theorem we prove a uniform Hausdorff 
dimension result for the level sets of B^'^ . 

Theorem 4.3 If 1/{HK) > d, then with probability one, 

dim^Z^ = 1 - HKd for all x G M"', (4.7) 

where dim^ denotes Hausdorff dimension. 



Proof It follows from Theorem 3.19 in Xiao (2005) that with probability one, 

dimjjZ^ = 1 - HKd for ah x € O, (4.8) 
where O is the random open set defined by 

0= y jxelR'^: L{x,[s,t]) >0}. 

s.tGQ; s<t 

Hence it only remains to show O = M.'^ a.s. For this purpose, we consider the stationary 
Gaussian process Y = {Y{t),t G M} defined by Y{t) = e~^^^B^'^ (e^), using the Lamperti 
transformation. 

Note that the component processes of Y are independent and, as shown in the proof of 
Proposition 12. H they are strongly locally (^-nondeterministic with (p{r) = r^^^. It follows 
from Theorem 3.14 in Xiao (2005) that Y has a jointly continuous local time LY{x,t), where 
{x,t) G M*^ X M. From the proof of Proposition 12. H it can be verified that Y satisfies the 
conditions of Theorem 2 in Monrad and Pitt (1987), it follows that almost surely for every 
y G M'^, there exists a finite interval J C M such that Lyiy, J) > 0. 

On the other hand, by using the occupation density formula (|4.ip . we can verify that the 
local times of B^'^ and Y are related by the following equation: for all 2; G M'^ and finite 
interval / = [a, b] C [0, 00), 

L{x,I)= I e^^-^^> LY{e-^^' X, ds). (4.9) 

J [log a, log b] 

Hence, there exists a.s. a finite interval / such that L(0, /) > 0. The continuity of L[x, I) 
implies the a.s. existence of 5 > such that L{y,I) > for all y G M'^ with \y\ < 6. 
Observe that the scaling property of B^'^ implies that for all constants c > 0, the scaled 
local time c''^^'^^'^^ L{x, ct) is a version of L(c ^^x, t). It follows that a.s. for every x G M'^, 
L{x,J) > for some finite interval J C [0, 00). □ 
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Since there is little knowledge on the explicit distribution of L(0, 1), it is of interest in 
estimating the tail probability P{L(0, 1) > x} as x — > oo. This problem has been considered 
by Kasahara et al. (1999) for certain fractional Brownian motion and by Xiao (2005) for a 
large class of Gaussian processes. Our next result is a consequence of Theorem 3.20 in Xiao 
(2005). 

Theorem 4.4 Let B^'^ = {B^'^ (t), t € M} be a bifractional Brownian motion in with 
indices H and K. If 1/{HK) > d, then for x > large enough, 

- log P{L(0, 1) > x} X x^^, (4.10) 

where a{x) x 6(x) means a(x)/6(x) is bounded from below and above for x large enough. 



Proof By Proposition 12.11 and Lemma 12.31 we see that the conditions of Theorem 3.20 in 
Xiao (2005) are satisfied. This proves (|4l(l . □ 

Let us also note that the existence of the jointly continuous version of the local time and 
the self-similarity allow us to prove the following renormalization result. The case d = 1 has 
been proved in Russo and Tudor (2006). 

Proposition 4.5 // 1/{HK) > d, then for any integrable function F : — > M, 

^HK-i I f(^b",k^^-^\ ^^(0, 1) as t^ oo, (4.11) 

J[o,t] 

where F = J^^ F{x) dx. 



Proof It holds that 

/ F{B''^''{u))du = t [ F{B'''''{tv))dv = t [ F (t^^5^'^(i;)) dz;. (4.12) 
Jlo,t] 7(0,1] J[o,i] 

By using the occupation density formula, we derive 

f F{B^'^{u))du = t [ F{t^^x)L{x,l)dx = t^-"^ [ F{y)L{yt-^^,l)dy. (4.13) 
Jlo,t] Jr'' Jm.'i 

Since y ^ L{y, 1) is almost surely continuous and bounded, the dominated convergence 
theorem implies that, as t — > oo, the last integral in (j4.13p tends to F L{0, 1) almost surely. 
This and iKWi yield (|iJT]) . □ 



4.2 Oscillation of bifractional Brownian motion 

The oscillations of certain classes of stochastic processes, especially Gaussian processes, in 
the measure space ([0, l],Ai), where Ai is the Lebesgue measure in M, have been studied, 
among others, by Wschebor (1992) and Azai's and Wschebor (1996). The following is an 
analogous result for bifractional Brownian motion. 
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Proposition 4.6 Let B^'^ be a bi-fBm in M. with indices H G (0, 1) and K € (0, 1]. For 

every t G [0, 1], let 

B^'^{t + e)-B^^^{t) 



Then the following statements hold: 
(i) For every integer k > 1, almost surely, 



[\z,it))''dt^E{p' 
Jo 



as e ^ 0, 



where p is a centered normal random variable with variance a"^ = 2^ ^ . 
(a) For every interval J C [0, 1], almost surely, for every a; € M 

Xi{t G J : Ze{t) <x}^ Xi{J)F{p <x) as e^O. 

Proof Let us denote 

y"'*= = / (ZJt))'' dt. 



Jo 



10 

It is sufficient to prove that 

Var(y^''=) < c{k) for some c{k) and (3 > 0. (4.14) 

Then the conclusions (i) and (ii) will follow as in Azai's and Wschebor (1996) by the means 
of a Borel-Cantelli argument. 

Note that ^ ^ 

Var(y^''=) = / [ Coy{Z,{u)'', Zsiu)'')dudv. 
Jo Jo 

We will make use of the fact that for a centered Gaussian vector (U, V), 



Cov(^7^ V'') = c{p,k)[Cov{U, V)Y [Var(C/)Var(y)] 

i<p<fc 



k—p 



Since the random variable has clearly bounded variance [cf . Lemma 12. 3| , it suffices to 
show that for every 1 < p < k, 



7' 

Jo 



[E{Ze{u)Ze{v))Y dudv < c,_g (4.15) 



We can write 



JO JO JO 



[E{Ze{u)Z,iv))Ydvdu = 2 / l^^_,^,-)[E{Z,iu)Ze{v))Y dvdu 

Jo Jo 

pi pu 

+ 2 / l^^_,>,)[E{Z,{u)Z,{v))Y dvdu 
Jo Jo 



■.= A + B. 
15 



Clearly A < ce, hence it suffices to bound the term B. Note that 

1 r r d^R 



E{Z,{u)Z,{v)) 



e 



2HK 



M— e J v—e 



dadb 



dbda. 



Since 



d'^R , , , 2HK 
-[a,b) 



dadb 



2^ 



2H , u2H\K-^ 2H-lu2H-l 



i2HK -l)\a-b\ 



2HK-2 



we have 



B < c{p, H, K) 



1 fu—e 



JO 



1 



-2HK 



+ c{p,H,K) [ I 
Jo Jo 

■.= Bi + B2. 



1 ru—e 



^a2H ^^2H^K-2 ^2H-1^2H-1 ^^^^ 



u rv 



dvdu 



1 



r2HK 



\a - U'^^^^'^dbda 



u—s J v—e 



dvdu 



The term B2 can be treated as in the fBm case [see Azais and Wschebor (1996), Proposition 
2.1] and we get B2 < ce^ for some constant /? > 0. Finally, since a?^^ + b'^^^ > a^^b^^ , 
we can write 

fl ru—e 

Bi<c{p,H,K) 



Jo 



1 



2HK 



u pv 



aHK-l^HK-l ^^^^ 



u—e J v—e 



dvdu 



c{p,H,K) 



1 ru-e ,^HK _ _ ^-^HKy ^ ^HK _ _ ^^HK^P 



Jo 



-HK 



rHK 



dvdu 



< c 



1 /^HK 



(u-e) 



HK\ P 



-HK 



dvdu 



A change of variable shows that Bi < ce^*-^ hk) ^ Combining the above yields (j4.15p . There- 
fore, we have proved (j4.14p . and the proposition. □ 



The above result can be extended to obtain the almost sure weak approximation of the 
occupation measure of the bi-fBm B^'^ by means of normalized number of crossing of B^'^ , 
where B^'^ represents the convolution of B^'^ with an approximation of the identity (t) = 

{^) with $ = l[-i,o]- If 

51 is a real function defined on an interval /, then the number of 

crossing of level u is 

Nu{9,l) = #{tel,g{t)=u}, 
where i^E denotes the cardinality of E. 

Proposition 4.7 Almost surely for every continuous function f and for every bounded in- 
terval I C M+, 

(-) e^-^^' / /(n)iV„(Sf'^,/)dn^ / f{u)L{u,I)du as e ^ 0. 
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Proof The arguments in Azai's and Wschebor (1996), Section 5, apply. Details are left to 
the reader. □ 



4.3 The multi-parameter case 

For any given vectors H = {Hi,...,Hn) G (0,1)^ and K = {Ki,...,Kn) G (0,1]^, an 
(iV, d)-bifractional Brownian sheet B^'^ = {B"'^{t),t G M^} is a centered Gaussian random 
field in R"' with i.i.d. components whose covariance functions are given by 



^ 1 r 



If ■ — Q A'^^j^o 



(4.16) 



It follows from (j4.16p that, similar to an [N, d)-fractional Brownian sheet [cf. Xiao and Zhang 
(2002), Ayache and Xiao (2005)], B^'^ is operator-self-similar. However, it does not have 
convenient stochastic integral representations which have played essential roles in the studies 
of fractional Brownian sheets. Nevertheless, we will prove that the sample path properties 
of B^'^ are very similar to those of fractional Brownian sheets, and we can describe the 
anisotropic properties of B^'^ in terms of the vectors H and K. 
We start with the following useful lemma. 

Lemma 4.8 For any e > 0, there exist positive and finite constants c^^ and c^g such that 
for all s,t e[e,l]^ , 



N 



.El 



Sj - tj\^"^^^ < E 



Bi (s) - B^ (t) 



N 



J ''J I ■ 



(4.17) 



and 



N N 
C4,r E - ^ detCov(sf '^(s), 5f '^(t)) < c,,3 \s, - tj 



(4.18) 



Here and in the sequel, detCov denotes determinant of the covariance matrix. 



Proof We will make use of the following easily verifiable fact: For any Gaussian random 
vector {Zi, Z2), 

detCov(Zi, Z2) = Var(Zi)Var(Z2|Zi), (4.19) 

where Var(2'i) and Var(Z2|Zi) denote the variance of Zi and the conditional variance of Z2, 
given Zi, respectively. 

By KW\i we see that for all s, t G [e, 1]^, 



detCov(5f'^(s),Bf'^(t) 



= E 
< E 



B^'^isf 
B^'^isf 



Var(i?f'^'(t)|i?f'^(s; 



E 



Bf'^(,s)-sr'"(t) 



H,K , 



(4.20) 
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Since Var (i?||^'^(s)) is bounded from above and below by positive and finite constants, it is 
sufficient to prove the upper bound in (j4.17p and the lower bound in (j4.18p . 

When N = 1, Lemma [231 Proposition O and imply that both ([CTT]) and 

hold. Next we show that, if the lemma holds for any B^'^ with at most n parameters, then 
it holds for B^'^ with n + 1 parameters. 

We verify the upper bound in (j4.17p first. For any s,t £ [e, 1]"^^, let s' = (si, . . . , s„, in+i)- 
Then we have 



E 



Bi (s) 



< 2E 



Bi [s) 



B^^^s') 



+ 2E 



(4.21) 



For the first term, we note that whenever si, . . . , s„ G [e, 1] are fixed, B^'^ is a (rescaled) 
bifractional Brownian motion in 5^=1 • Hence Lemma 12.31 implies the first term in the right- 
hand side of (j4.2ip is bounded by c|t„ — where the constant c is independent 
of si,... ,s„ G [£,1]. On the other hand, when tn+i ^ [^i^ is fixed, B^'^ is a (rescaled) 
(iV, (i)-bifractional Brownian sheet. Hence the induction hypothesis implies the second term 
in the right-hand side of (I4.21j) is bounded by c Yl^=i l*i ~ SjP^^^^ - This and (I4.2ip together 
prove the upper bound in (j4.17p . 

Suppose the lower bound in ()4.18p holds for any B^'^ with at most n parameters. For 
= n + 1, we write detCov(i?^'^(s), B^'^ {t)) as 



n+l 



n+l 



2HjKj 2HjKj 



n+l 

n 



1 



22K, [_V3 



j=2 



^HjKj 2HjKj \ 2HiKi .2HiKi 



'1 



1 



22Ki [V I 



(t2H, ^ ^^2H,^K, 



\ti - Si 



\2HiKi 



+ 



1 



22Ki 

n+l 



(tf. +sf')*"'-|(i-s,p».'<-> 



2HjKj 2HjKj 



j=2 



n+l 

n 

j=2 



22^^ 



tp + spy^-\tj 



2Hi \ Ki 



n+l 



bj ij I 



(4.22) 



where the last inequality follows from the induction hypothesis. This proves the lower bound 
in (liT8|) . □ 

Applying Lemma 14.81 we can prove that many results in Xiao and Zhang (2002), Ayache 
and Xiao (2005) on sample path properties of fractional Brownian sheet also hold for B^'^ 
as well. Theorem 14.91 is concerned with the existence of local times of B^'^ . 



18 



Theorem 4.9 Let B^'^ = {S^'^(t), t S M^j he an {N, d)-bifractional Brownian sheet with 
parameters H G (0,1)^ and K G (0,1]^. If d < Y.f=nrK- ^^^n for any N -dimensional 
closed interval I C (0,oo)''^, B^'^ has a local time L{x,I), x G M.'^. Moreover, the local time 
admits the following L'^ -representation 

L{x,I) = {27r)-^ [ e-^'^y^''^ f e'^y^^"-''^'^Usdy, x € M'^. (4.23) 

Remark 4.10 Although the existence of local times can also be proved by using the Malliavin 
calculus [see Proposition 14.151 below] , we prefer to provide a Fourier analytic proof because: 
1) we can compare in this way the two methods and 2) the above theorem gives in addition 
the representation (I4.23[) . 

Proof Without loss of generality, we may assume that I = [e, 1]^ where e > 0. Let Atv be 
the Lebesgue measure on /. We denote by /x the image measure of A at under the mapping 
1 1-^ B^'^[t). Then the Fourier transform of /i is 

/i(^) = ^e^(€.^"'"W)dt. (4.24) 

It follows from Fubini's theorem and (j4.17p that 

e[ \jl{0\'d^=[ [ [ E(e^«'^"'"(^)-^"'"W))dCd.dt 

jR'i J I J I JRd ^ ^ 

^'//[E(Bf'^(.)-V'^(t))Y/^^'^' (4.25) 
< c [ [ TTTT dsdt. 

The same argument in Xiao and Zhang (2002, p. 214) shows that the last integral is finite 
whenever d < X^jLi jtk'- Hence, in this case, fi G L'^{W^) a.s. and Theorem 14. 91 follows from 
the Plancherel theorem. □ 



Remark 4.11 Recently, Ayache, Wu and Xiao (2006) have shown that fractional Brownian 
sheets have jointly continuous local times based on the "sectorial local nondeterminism" . 
It would be interesting to prove that B^'^ is sectorially locally nondeterministic and to 
establish joint continuity and sharp Holder conditions for the local times of B^'^ . 

Now we consider the Hausdorff and packing dimensions of the image, graph and level set 
of B^'^ . In order to state our theorems conveniently, we assume 

< HiKi <...< HnKn < 1. (4.26) 

We denote packing dimension by dimp; see Falconer (1990) for its definition and proper- 
ties. The following theorems can be proved by using Lemma 14.81 and the same arguments as 
in Ayache and Xiao (2005, Section 3). We leave the details to the interested reader. 
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Theorem 4.12 With probability 1, 

N 



dim,5^'^([0, 1]^) = dim,S^.^([0, 1]^) = min U E 1^1 (4-27) 



and 



dim„GrB^-^([0, 1]^) = dimpGrS^-^([0, 1]^) 



^ f EjLl HjKj ^f=l H^K, ^ 

\T.U^ + N-k + {l-H,K,)d if E,tJH^<ci<E,ti^, 

(4.28) 

where E?=i Hjjq ■= 0- 

Theorem 4.13 Let = {t e (0, oo)^ : B^^^{t) = x} be the level set of B^'^ . The 
following statements hold: 

(i) //E^Li w ^ ^' ^^^^ /^^^ every x ElW^ we have = $ a.s. 
(a) If EjLi > ^; then for any x eW^ and < e < 1, with positive probability 



dim^ (L^ n [e, If) = dinip (L, n [e, if) 

k 

= mm 



l^^+iV-A: - Hkd, 1 < A; < 

k -w-T k 1 ^ k ^ 



(4.29) 



4.4 A Malliavin calculus approach 

Using the Malliavin calculus approach, we can study the local times of more general bifr ac- 
tional Brownian sheets. Consider the (A'' x d)-matrices 

H = (Hi,...,Hd) and K = {Ki, . . . ,Kd), 

where for any i = 1, . . . ,d 

Hi = {Hi,u Hi^N) and Ki = {Ki^i, Ki^N) 

with Hi^j G (0, 1) and Kij G (0, 1] for every i = 1, . . . , c? and j = 1, . . . ,N. 
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We will say that the Gaussian field B^'^ is an {N, d)-bifractional Brownian sheet with 
indices H and K if 

and for every i = l,...,d, the process {B^^{t), t G R^} is centered and has covariance 
function 

N 

i=i 

As in subsection 4.1, the local time L{x, t) { t G and x € W^) of B^'^ is defined as 
the density of the occupation measure fXt, defined by 

/xt(A) = / ds, A e BiR'^). 

J[0,t] 

Formally, we can write 

L{x,t)= [ 5^{B'^^{s))ds, 

J[Q,t\ 

where Sx denotes the Dirac function and Sx {Bg ' ) is therefore a distribution in the Watanabe 
sense (see Watanabe (1984)). 

We need some notation. For a; € M, let Pa{x) be the centered Gaussian kernel with 
variance a > 0. Consider also the Gaussian kernel on M'^ given by 

d 

Pt{x) =Wp„{xi), x = {xi,...,Xd) eR'^. 
1=1 

Denote by H„(x) the n-th Hermite polynomial defined by Ho(a;) = 1 and for n > 1, 

H„(x) = exp (-) — exp (--) , x G R. 

We will make use of the following technical lemma. 

Lemma 4.14 For any H G (0, 1) and K G (0, 1], let us define the function 

^ , , i?^'^(l,z) , , 

Qh,k{z)= ^ ^G(0,1] 

and QH,KiO) = 0. Then the function Qh,k takes values in [0,1], Q_f/,x(l) = 1 (^nd it is 
strictly increasing. Moreover, there exists 6 > such that for all z & {1 — 6, 1), 

{QhM^T < exp (-c(<5, H, K)n{l - zf") . (4.30) 
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Proof Clearly, the Cauchy-Schwarz inequality implies < Qh,k{z) < 1- Let us prove 
that the function Qh,k is strictly increasing. By computing the derivative Q'jj x{z) and 
multiplying this by z^^~^^, we observe that this is equivalent to show 

{l-zf"^-\l + z)-{l + z'^")^-\l-z^^)>0 forall zG(0,l). (4.31) 

If HK < t;, since (1 + < 1 + z, the left side in (ITOI) can be minorized by 

(1 + ((1 - - 1 + z^") and this is positive since (1 - > 1. 

If HK > ^, we note that 

(1 _ zftiK^^l + z) + (1 + > (1 - Z)(l + Z) + (1 + Z2^^)^-Iz2 

> (1 + ^2^^)^-1(1 _ z^) + (1 + Z^"f-'Z^ > (1 + Z^'^f-K 

and this implies (j4.3ip . Concerning the inequality (j4.30p . we note that 

Qh,k{zT = exp {n\ogQH,K{z)) > exp (-n(l - Qh,k{z))) ■ 
Now by Taylor's formula 

(1 + z'^'fz-'''' < 2^ + c{H, K, d){l - zf 

and therefore 

Qh,k{^) < 1 + ciH,K,6)il - zf - ^(1 - zf"^ 

< 1 + C{H, K, 5){l - z)2^^^j2-2HK - -L(l _ ^)2Hi^ . 

The conclusion follows as in the proof of Lemma 2 in Eddahbi et al. (2005), since 

1 - QhM^) > - - 6)) 

for any z E (1 — 5, 1) with 5 close to zero and with c{H, K, 6) tending to zero as 5 — > 0. □ 



The following proposition gives a chaotic expansion of the local time of the {N, d)- 
bifractional Brownian sheet. The stochastic integral In{h) appeared below is the multiple 
Wiener-Ito integral of order n of the function h of nN variables with respect to an (N, 1) 
bifractional Brownian motion with parameters H = {Hi, . . . , Hj^) and K = {Ki, . . . , K]\j) . 
Recall that such integrals can be constructed in general on a Gaussian space [see, for example. 
Major (1981), or Nualart (1995)]. We will only need the following isometry formula: 

AT 

E (/„(l|p/„(]lfo™ )) = i?^'^(t, sT\n=m) = n (^''"'"^ ' '3)T \n=m) (4-32) 

i=i 

for all s,t G 
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Proposition 4.15 For any x G and t G (0, oo)^, the local times L{x, t) admits the 
following chaotic expansion 

i(-.*)= E L,n^fil^«".(f)^".(v.i(-)-)rf., (4.33) 

where s = si---sn and s^^^^ = Y\jLi s^""'^^''^ ■ The integrals denotes the multiple ltd 
stochastic integrals with respect to the independent N -parameter bifractional Brownian motion 

Moreover, if X^^i h*k* ^ '^^^''"^ H* = max{Hij : i = 1, . . . , d} and K* = m.ax{Kij : 
i = 1, . . . , d}, then L{x, t) is a random variable in L'^{Q). 



Proof The chaotic expression (|4.33p can be obtained similarly as in Eddahbi et al. (2005) 
or Russo and Tudor (2006). It is based on the approximation of the Dirac delta function 
by Gaussian kernels with variance converging to zero. Let us evaluate the L'^{Q) norm of 
L(x, t). By the independence of components and the isometry of multiple stochastic integrals, 
we obtain 



where 




(4.34) 



By Propositions 3 and 6 in Imkeller et al. (1995) [see also Lemma 11 in Eddahbi et al. 
(1995)], we have the bound 

/3n, {u)Pn, {v) < c, J- (4.35) 

for any (3 G [j, |). Using the inequality (|4.35p . we derive from (j4.34p that \\L{x, t)||| is at 
most 

^E E (n— T^)//"/ ,^"nn^7^^ 

^>Om + -+na=m^i=l{niVl) 6 / J[0,t] JlO,u] j=i {^3^3 ) 

= ^ E E (n ; — W) n t ^^-^^^^ /' (UQ^^..^^A^r)'^z (4.36) 

m>Oni+-+na=m^i=l [riiy 1) 6 / ^.^^^0 ^0 V / 

"■^....^E E (n— ^)n/'(n<3"...-...(-)"0*- 

m>OniH hna=m ^i=l [niV L) 6 / j^^JU \ -^^ / 
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where we used the change of variables uj = uj and vj = zjuj. Using the above lemma and 
as in the proof of Lemma 2 in Eddahbi et al. (2005), we can prove the bound 

I' (nQ/^..,i^..(^)"')rf^ < C4.n (4.37) 
Here = c^^^(i7, K) depends on H,K. Finally, (14.370 implies that 

N IX . d 



i(^.*)iii<c.,„E(n'»' ) E (n--^ 

in.; V L] 6 



m>l ^j=l ^ ni + ---+nd=m "^4=1 (^jVl) 6 



(4.38) 



m 



m>l 

where and jg depend on i/, K and t only. The last series in (|4.38p converges if 

f;^>rffl-^V (4.39) 

To conclude, observe that by choosing f3 close to |, XljLi jj-*^* ^ ^ implies the required 
condition (|i39D . □ 



We recall that a random variable -F = Inifn) belongs to the Watanabe space B"'^ if 

ll^ll',2 :=5](l+m)-||/„(/„)||2<oo. 

n>0 

Corollary 4.16 For any t G (0, oo)^ and x E , the local time L{x, t) of the {N,d)- 
bifractional Brownian sheet B^'^ belongs to the Watanabe space B"'^ for every < a < 

SpN 1 _ d 

2^j=l 2H*K* 2 • 

Proof This is a consequence of the proof of Proposition 14.151 Using the computation 
contained there, we obtain for any /3 S [5, ^), 

\\L{x,t)\\l,<c,MH,K,d,t)Y,il + mrm'^ « )" -^.^1 ^ 

m>l 

which is convergent if a < Ylf=i 2H*k* ~ ~ ^^6^ ) — 1 — Ylf=i 2H*k* • Choosing /3 close 

■' j 3 ■'33 

to i, we get the conclusion. □ 
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